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Abstract 

The properties of compact stars made of massive bosons with a repulsive selfinteraction mediated 
by vector mesons are studied within the mean-field approximation and general relativity. We 
demonstrate that there exists a scaling property for the mass-radius curve for arbitrary boson 
masses and interaction strengths which results in an universal mass-radius relation. The radius 
remains nearly constant for a wide range of compact star masses. The maximum stable mass and 
radius of boson stars are determined by the interaction strength and scale with the Landau mass 
and radius. Both, the maximum mass and the corresponding radius increase linearly with the 
interaction strength so that they can be radically different compared to the other families of boson 
stars where interactions are ignored. 

PACS numbers: 05.30.Jp Boson systems, 04.40.Dg Relativistic stars, 04.40.-b Self-gravitating systems, 
05.70.Ce Thermodynamic functions and equations of state, 95.35.+d Dark matter 
Keywords: compact stars, boson stars, interacting bosons, dark matter 
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I. INTRODUCTION 



White dwarfs, neutron and quark stars, collectively dubbed compact stars, are the final 
result of stellar evolution. White dwarfs are stabilized by the Fermi degeneracy pressure. 
There exists an upper limit for the mass of a white dwarf, the Chandrasekhar mass, which 

n 

is about 1.4 times the mass of the sun [l|. Beyond this limit the white dwarf is unstable 
against gravitational collapse. Neutron stars are stable mainly due to the repulsive nature of 
the interactions between nucleons. Therefore, the precise value of the maximum mass for a 
neutron star is less certain, but is presumably close to the predictions based on the Landau 

li n 

consideration [2|. As shown in Ref. [3|, Landau's argument can be extended to a general 
compact star made of fermions, a fermion star, with arbitrary fermion mass and interaction 
strength. 

In the following we are studying compact stars made of bosons. Unlike fermion stars, 
boson stars have no observational evidence, yet. Besides this the existence of any stable 
scalar particle has never been experimentally verified. Wheeler ^| introduced a gravitational 
electromagnetic entity called a geon. The gravitational attraction of its own field energy 
confines the geon in a certain region. Later Kaup [5] has solved the Klein-Gordon Einstein 
equations for scalar fields and found a new class of solutions for gravitating objects. These 
boson stars are stable with respect to spherically-symmetric gravitational collapse. Ruffini 
and Bonazzola [fj] demonstrated that boson stars describe a family of self-gravitating scalar 
field configurations within general relativity. In ref. [?| Takasugi and Yoshimura calculated 
boson stars within an approach similar to the one conventionally adopted for neutron stars by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equation j^, [h]] with a separate equation 
of state describing the properties of matter. Boson stars with selfinter act ions have been also 



considered^ in particular as candidates for dark matter ll|] . For reviews on boson stars we 



refer to 



In the following, we consider boson stars as localized, gravitationally bound objects made 
of self-interacting bosons at zero temperature. The interactions between the bosons is de- 
scribed by vector meson exchange in the relativistic mean-field approximation. The resulting 
equation of state is used as input to solve the TOV equation for boson stars, similar to the 
approach of ref. \j] but with an equation of state based on a field-theoretical approach. 
We demonstrate that there are scaling relations for the mass-radius curve. In particular, 
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we show that the maximum mass is controlled by the interaction strength and the Landau 
mass, not by the boson mass. We compare our results to previous works and to the case of 
fermion stars with self-interactions. 



II. SCALING RELATIONS FOR COMPACT STARS 

We assume a spherically-symmetric and static configuration where the energy-momentum 
tensor is that of a perfect fluid at rest. Then the star structure can be obtained by solving 
the Tolman-Oppenheimer-Volkoff (TOV) equations, which can be conveniently written as 

with 

_ _ = 4vrr 2 p . (2) 
dr 

Additionally, one needs an equation of state, p(p), which describes the microscopic properties 
of the stellar matter. These coupled differential equations for the pressure p(r) and the 
mass profile M(r) are integrated from r = with some central value po to a point where 
the pressure vanishes p(R) = which defines the radius R and the total mass M(R) of the 
boson star. 

The TOV equations have a similar scaling behavior as the one for Newtonian hydrostatic 
equilibrium (see e.g. Q]). Just the gravitational constant G and the boson mass m& are 
used to rewrite the TOV equations together with the equation of state in dimensionless 
form. Solving the dimensionless TOV equation for a certain class of equations of state 
allows for deriving general solutions by just rescaling the results by appropriate dimensionfull 
parameters. 

The first relativistic correction factor in eq. ([I]), i.e. (1 +p/p), can be scaled by choosing 
p' = Pi Po and p' = pj p - Here p Q is a common factor with dimension of mass to the fourth 
power. For the other two relativistic factors in eq. (CQ) we introduce the dimensionless mass 
M' = M/a and the dimensionless radius r' = r/b. For a dimensionless expression one has 
to set 



= ! M.I 1(1 

MU 



1 and THl = l ( 3 ) 
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that leads to the following relations 



Ml and b= ML. (4) 



where Mp = ^Jhc/G is the Planck mass. Note that the Newtonian terms do not give 
any additional constraint. Choosing p Q = mf, the rescaling factors are a = Mpjm\ and 
b = Mp/mf, which coincide with the expressions of the maximum mass and the radius for 
compact stars introduced by Landau 

M L = ^ and R L = ^ (5) 

We note in passing that for the MIT bag equation of state these scaling factors are p Q = B 
so that the maximum mass and the corresponding radius scale as 5 _1//2 , a well known result 
for quark stars [is! ]. 

III. MESON EXCHANGE MODEL FOR INTERACTING BOSONS 

We describe the interactions between scalar bosons by the exchange of vector mesons. 
For a scalar field and a vector field the Lagrangian reads 

C = Vl4TVy - m b y - + \mlV^ (6) 

with V^v = d^V v — d u Vfj,. The boson field is coupled to the vector field by a minimal coupling 
scheme 

1)^ = 9^ + ig^Vn . (7) 

where g v ^ is the (p-V coupling strength. Note, that the vector field has a quadratic coupling 
term to the scalar field in the Lagrangian which ensures that the vector field is coupled to 
a conserved current (see below). We treat the vector field as a classical field. In static bulk 
matter the spatial components of the vector field vanish and the equation of motion for the 
scalar field reads 

'V*V + m b 2 ] 4>(x) = . (8) 

In the mean-field approximation, after expanding into plane waves, we obtain for the lowest 
energy mode k = 0: 

Utf, = m b + g v4> V Q (9) 
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Note that the vector interaction between the scalar particles is repulsive which ensures the 
overall stability of selfinteracting boson matter. The vector field is determined from the 
equation 

m 2 v V = 2g V(j> (u^ - g v<t> V ) <fr*<p = 2g v4> m b (f)*(f) (10) 

where we have used the dispersion relation for the field eq. (Q. The conserved current for 
the scalar field can be obtained from the Lagrangian (jSJ) 

. / dC dC \ 

= fid^ - {id„<t>*)<t> + 2g V(t y^*^ . (11) 

The number density of bosons 

n b = J = 2 (uj, - g v <j>Vo) <f)*<f) = 2m b 4>*4> (12) 

is just the source term for the vector field. The total energy density of the boson matter 
can be determined from the energy-momentum tensor 

1 9 2 

p = 2m\<\)*(\) + -m 2 v Vo = m b n b + 7^ n & ( 13 ) 

where the equation of motion for the vector field has been used. The pressure is given just 
by the vector field contribution 

H m ^if> < 14 > 

Note, that these expressions are thermodynamically consistent as can be checked by using 
the thermodynamic relation 

D , ^(pM) (15) 

P ~ Ub dn b (15j 
The form of the interaction is actually similar to the one used for interacting fermions and 
the corresponding Fermi stars in [3]. For the rescaled TOV equations we introduce the 
dimensionless interaction parameter y = m b /mj, where m/ = \/2m»/^, and set p' = pjm\ 
and p' = p/m b . The equation of state for interacting boson matter can be summarized to 
be of the simple form 

p' = y 2n 'b and P = n 'b + y 2n 'b ( 16 ) 

with the dimensionless number density n' b = n b jvn^. It is possible to represent the equation 
of state in a polytropic form p = p 1 for certain limits. For low densities, one approaches 
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FIG. 1: Double logarithmic plot of the dimensionless pressure versus the dimensionless energy 
density for different interaction strengths y. 

p oc p 2 , a polytrope with 7 = 2. For high densities, one has an equation of state of the form 



P 



p with 7 = 1, which is the stiffest possible equation of state first discussed by Zeldovich 



161 ] . The switch between those two limiting cases is controlled by the interaction strength 
y. The larger y, the lower is the energy density to approach the causal limit p = p. 

Fig. [1] depicts the dimensionless pressure versus the dimensionless energy density. The 
values of the interaction strength y are chosen between 10~ 2 to 10 7 . There are two different 
slopes for small and large values of p' corresponding to the above mentioned limits. The 
point where the slope changes shifts to lower densities with increasing interaction strength 
y, but for low densities the slope is 7 = 2. At high p' all curves merge to the limiting curve 
with a slope of 7 = 1. 



IV. SCALING RELATION FOR BOSON STARS WITH SELFINTERACTION 

We use the dimensionless equation of state to solve the dimensionless TOV equations. 
The equation of state depends only on the interaction strength y. Fig. [2] shows the double 
logarithmic plot of the dimensionless mass M' versus the dimensionless radius R' for different 
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FIG. 2: The dimensionless mass M' is plotted versus the dimensionless radius R' for different 
interaction strengths y ranging from 10~ 3 to 10 5 . Note that each mass-radius curve terminates in 
a spiral at the left end, which is not visible in the double logarithmic plot, but can be seen in the 
linear plot of Fig. [H 

interaction strengths y ranging from 10~ 3 to 10 5 . One observes that each mass-radius curve 
contains a constant radius part over a wide range of masses. Also, the curves are very similar 
and seem to be just shifted to larger masses and radii with increasing interaction strength. 

This interesting behavior can be explained by considering the equation of state described 
by a polytrope p ~ p 7 with 7 = 2. In general, the solution to the Lane-Emden equation, 
see e.g. [l7j, results in a mass-radius relation of the form M' oc pf 1-4 ^ 2 and R' oc p^ 1 ^ 2 ^ 2 , 
where p c is the central energy density. Hence, at low densities and large radius R 3> 2GM, 
where effects from general relativity can be ignored, the mass of the star increases linearly 
with p c while the radius R' remains constant for 7 = 2 that explains the peculiar form of 
the mass-radius curves. 

There exist another interesting feature of the mass-radius curves which reflects the scaling 
properties of the equation of state and the TOV equations. To illustrate this we plot in Fig. [3] 
the dimensionless maximum mass as a function of the interaction strength y. It is 
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FIG. 3: The dimensionless maximum mass M max is plotted as a function of interaction strength y. 

interesting to see that the maximum mass M max scales linearly with the interaction strength 
V- 

Therefore, one can conclude that by proper rescaling all mass-radius curves can be reduced 
to one universal mass-radius curve. Indeed, dividing the dimensionless mass M' and the 
corresponding radius R' by the interaction strength y results in an unique mass-radius 
relation as depicted in Fig. HI This graph looks rather similar to the mass-radius curve 
of a strongly interacting fermion star [3|. There the maximum mass is constant for weak 
interactions (y< 1) and increases linearly in y for strong interactions (y 3> 1). Note that 
the part of the curve to the left of the maximum mass represents unstable configurations, 
only the star configurations at the maximum and to the right of it can exist. The existence 
of a maximum mass is determined by the change of the equation of state from a polytrope 
with 7 = 2 to one with 7 = 1. The latter value is lower than the critical value 7 C = 4/3 for 
stable compact stars (effects of general relativity will even increase this value slightly). 

Fig. [5] shows the two limiting radii for interacting boson stars as a function of the interac- 
tion strength y. Here, -R max and R' min denote the maximum and minimum radius for boson 
stars, respectively. R' min stands for radius corresponding to the maximum mass configura- 
tion, while -R max is the radius of stars with masses much smaller than the maximum mass. 
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FIG. 4: The dimensionless mass M' is plotted versus the dimensionless radius R' dividing both by 
the interaction strength y. The result is an universal mass-radius relation. 





FIG. 5: Plot of the two limiting radii R' max and R' min for boson stars as a function of interaction 
strength y. 
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FIG. 6: The variation of the normalized density p'(r)/p'(0) with the normalized radius r'/R' for 
different interaction strengths showing that there is only one universal curve. 

Both, the maximum and the minimum radius vary linearly with the interaction strength y 
and the difference between the two radii is rather small, by a factor of about 0.61 independent 
of the interaction strength. 

Fig. [6] shows the normalized density profile of p'{r)/p'(0) over the normalized radius r'/R' 
for different interaction strengths calculated for the maximum mass configuration. Again, 
there appears an universal curve independent of the interaction strength and the mass of 
the boson. The rate of the decrease of the density with the radius is then the same for all 
interaction strengths y. The density profile shows a small plateau in the core region up to a 
radius of about r ~ O.li? followed by a nearly linear decrease up to the surface of the boson 
star. 

The scaling behavior observed for interacting boson stars follows in straightforward way 
from our discussion above on the scaling features of the TOV equations. The maximum mass 
configuration is determined by the equation of state at high central densities, where p = p. 
The pressure and the energy density, p and p, depend then on the interaction strength as y 2 . 
Rescaling both, the energy density and the pressure, by the factor mf/y 2 gives the modified 
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Landau mass and Landau radius 

Ml = y^i and R b L = y ■ ^ (17) 

m^ mj; 

for compact stars with interacting bosons. The maximum mass and the corresponding radius 
have to increase linearly with the interaction strength y. This is in complete agreement with 
the results of Ref. 



llj where the selfmteraction term of the form A0 4 was used to describe 
the interactions between bosons. Also for this type of interactions, the maximum mass 
and the corresponding radius were found to scale with the interaction strength and the 
Landau mass as M max oc \ 1/2 M L . The relation for the maximum mass is compatible with 
our findings by realizing that the dimensionless coupling constant A can be associated with 
our interaction strength y 2 . In addition to the case of scalar selfmteraction, we find that for 
vector interactions the scaling property is even more general as the whole mass-radius curve 
can be described by an universal curve when using the modified Landau mass and radius. 



We note, that compact stars made of fermions with vector interactions |3] reveal the same 
scaling feature of the mass-radius curve for large interaction strengths. 

The maximum mass of boson stars as obtained from numerical calculations is 

M 3 

M max ^0.16%-— J (18) 
m 2 

and the two limiting radii of boson stars are given by the expressions 

iU « 1.2523/ • — f and R min « 0.763?/ • — £ (19) 
mi mi 

The above relations can be used to calculate the maximum mass and the maximum and 
minimum radii of boson stars for arbitrary interaction strength y and boson masses mj. 
The values for the Landau mass and radius are Ml = 1.632M Q and Rl = 2.410 km, 
respectively, for a boson mass of m& = 1 GeV. One recovers the same scaling relations as for 
the noninteracting case, see e.g. Ref. by setting the interaction scale equal to the Planck 
mass, mi = Mp, 

M 2 1 
M max oc — - and R min oc — , (20) 
m b m b 

which are orders of magnitude smaller than for the case of realistic interactions. We note that 
our numerical prefactors are different from the ones of Takasugi and Yoshimura [jj while the 
scaling with the boson mass is the same. These authors adopt a different equation of state, 
where the pressure has the form as for degenerate configurations, e.g. in the low-density 
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y = 1 


R (km) 


10- 4 
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(Landau limit) 


M max (M ) 


10~ 5 


10~ 3 
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10 3 


10 6 


10 11 


mi = 100 GeV 


R (km) 


10- 4 


10" 2 


10 
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(weak scale) 


M max (M ) 


10- 2 


1 


10 3 


10 6 


10 9 


10 14 


mi = 100 MeV 


R (km) 


0.1 


10 


10 4 


10 7 


10 10 


10 15 


(QCD scale) 



TABLE I: Order of magnitude scales of the maximum mass and the characteristic radius of 
compact stars made of different boson masses and interaction strengths. The first set corresponds 
to the free case by setting y = m^/Mp (the boson stars are just bound by gravity), the second 
gives the Landau mass and radius by setting y = 1, which holds for boson and fermion stars. The 
third and last set lists the values for interaction mass scales of the standard model weak and strong 
interactions, i.e. of 100 GeV and 100 MeV, respectively. 

limit they recover that p oc p 5 / 3 . In our case the pressure is determined by interactions only 
and is proportional to the density squared. 

Table d gives the maximum mass and the corresponding radius for four different cases 
of the interaction parameter. By setting mj = Mp or y = m^/Mp one recovers the case 
for ordinary boson stars with free bosons (see above) which are just gravitationaliy bound. 
The case y = 1 gives the Landau mass and radius of compact stars, which is nearly the 
same for boson stars and fermion stars when including interactions. Finally, we consider the 
case of interactions mediated by the weak interaction scale of about mi = 100 GeV and the 
QCD scale of about mi = 100 MeV. For the boson masses we choose the range from the 
electroweak scale to a typical mass of axions, ~ 10 -5 eV. We want to emphasize the following 
features of these calculations. For the free case mi = Mp one only reaches astrophysically 
interesting scales for boson masses of less than 10 -5 eV. Maximum masses close to the ones 
for neutron stars, M ~ 1M , can be reached by boson stars for boson masses of around 
1 GeV (Landau case), 1 MeV for bosons with weak interactions, and 1 GeV for strong 
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interactions. The mass range of observed supermassive black holes, M = 10 6 to 1O 9 M , is 
found for boson masses between 1 keV and 1 MeV for the Landau case, 1 eV and below for 
weak interactions, and 1 eV to 1 keV for strong interactions. It is observed that the inclusion 
of interactions results in a wide range of possible masses and radii for boson stars, covering 
scales as small as a fraction of a solar mass and below a kilometer to scales of supermassive 
black holes. It is interesting to note, that a boson with a mass of 100 GeV and with QCD- 
type interaction strengths gives star configurations with masses and radii as a neutron star. 
Surprisingly, for a boson star made of axions (mj ~ 10~ 5 eV) and selfinteractions on the 
scale of 10 12 GeV, one obtains a mass of about 3OM with a radius of 200 km, i.e. the 
mass of compact objects found in binary systems which are attributed to light black holes. 
These values are orders of magnitude different compared to the case of boson stars with 
noninteracting axions, see also ref. [7| and table [B 



V. SUMMARY 



We have constructed an equation of state for a system of massive bosons interacting by 
the exchange of vector mesons. By solving the TOV equations for such boson stars, we have 
demonstrated that there exists a universal mass-radius curve independent of the boson mass 
and the interaction strength. The maximum mass and the corresponding radius of boson 
stars are scaled with the Landau mass and Landau radius times the interaction strength. 
For masses much smaller than the maximum mass, the radius stays constant and is only 
slightly larger than the one for the maximum mass configuration. 

The maximum mass and the corresponding radius can be computed with the simple 
formulae M max = 0.16% ■ Mp/mf and -R m i n = 0.763?/ ■ Mp/ml for any given boson mass 
rrib and interaction strength y. The possible masses and radii for boson stars can therefore 
cover a wide range and can be similar to the ones found for astrophysical compact objects, 
be it neutron stars or black hole candidates. For example, for a boson with QCD-type 
interaction strength and a boson mass of 100 GeV the maximum mass is M max ~ O.3M 
with a radius of about 2 km. A boson with a typical axion-like mass of 10 -5 eV and an 
interaction scale of 10 12 GeV will give a maximum mass of the boson star of 3OM with a 
radius of 200 km. The compactness of boson stars for the maximum mass configuration is 
about R/(2GM) w 2.3 which is close to the value found for fermion stars R/(2GM) « 2.4 
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in Ref. |3j. It is interesting that these values are below the radius of the innermost stable 
circular orbit of nonrotating black holes R/(2GM) = 3. 

Finally, we mention that the full problem addressed here involves solving the Einstein 
equations with the coupled system of Klein-Gordon and Proca equations which we leave to 
address as an interesting extension for future work. 
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